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1. Introduction 

In the context of quantum gravity, it has been often advocated the possibihty that higher-derivative 
corrections should be added to the Einstein- Hilbert (EH) lagrangian of general relativity, in order to obtain 
a better behaved theory. Among the various possibilities, an especially interesting class is given by the 
so-called Lovelock Lagrangians. 

These were introduced in [1] as the only possible generalization to higher dimensions of the EH lagrangian 
that gives rise to field equations which are second order in the metric, linear in the second derivatives, and 
divergence free. 

Each term Li of the Lovelock Lagrangian consists in a combination of contractions of a given power of the 
Riemann tensor, the lowest order terms Lq and Li being the cosmological constant and the EH Lagrangian, 
respectively. The i-th term is defined in dimensions d > 2i, and in d = 2i reduces to a total derivative, that 
takes the form of the Gauss-Bonnet topological invariant. For this reason, the Lovelock Lagrangians are 
sometimes called Gauss-Bonnet (GB) lagrangians. In the following, for simplicity of notation, we reserve 
this name for the quadratic term 

which gives the lowest order correction to the EH lagrangian and hence is the most common in applications. 
This quadratic combination was originally introduced by Lanczos [2] . 

Prom the topological origin of the Lovelock terms it also follows the important property that they do 
not introduce any new propagating degrees of freedom in the spectrum [3]. In particular, L2 is the only 
quadratic combination of curvature terms that does not give rise to tachyons or ghosts. In four dimensions 
it reduces to a total derivative, and hence does not contribute to the field equations unless nonminimally 
coupled to a scalar field. 

Because of their properties, the Lovelock Lagrangians were widely investigated in the context of Kaluza- 
Klein theories [4]. Later it was realized that they intervene also in the low-energy limit of string theories, 
where the GB term, nonminimally coupled to the dilaton, gives a contribution to the effective action [3,5]. 

The presence of higher-order terms in the action may smooth the behaviour of the solutions in the 
regime of strong gravity, as for example near the initial singularity of cosmology, and can induce spontaneous 
compactification. Cosmological solutions of EH-GB models were therefore investigated in [6]. More recently, 
a renewed interest in Gauss-Bonnet cosmology has arisen in the context of the braneworld scenario [7]. 

Lovelock Lagrangians are also of interest for black hole physics. In fact, in addition to the modification 
of the short-distance behaviour of the solutions near the singularities, also global properties may be altered 
if one allows for the possibility of nonminimal coupling of scalars as required by the string effective action, 
since the standard no-hair theorems [8] can be circumvented in that case. 

Exact black hole solutions of the higher- dimensional EH-GB model in the absence of scalar fields were 
obtained in [9], and exhibit either Minkowski or (anti-)de Sitter asymptotics. The more difficult case of 
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the EH-GB theory coupled to a scalar field in four dimensions was originally studied perturbatively in [10], 
using the methods introduced in [11] for higher dimensions, and then numerically in [12]. It turned out 
that regular asymptotically flat black hole solutions exist with nontrivial scalar fleld. The scalar charge is 
however a function of the mass (and charge) of the black hole'f. 

In these papers the existence of black hole solutions that are not asymptotically flat was not investigated. 
However, it is known that in D > 4 the EH-GB theory with no scalar possesses (anti-)de Sitter solutions 
[3-4]. Although in [14] it was argued that no asymptotically (anti-)de Sitter solutions can be present if a 
nonminimally coupled scalar fleld is added to the model, it is interesting to seek for more general solutions 
of the EH-GB system, with arbitrary asymptotic behaviour. 

A powerful technique for investigating this topic is the study of the phase space of the solutions of 
the fleld equations. This method has been used in several models of Einstein-Maxwell theory nonminimally 
coupled to scalars [15]. As mentioned above, when a GB term is added to the action, the fleld equations 
are still second order, and linear in the second derivatives, but no longer quadratic in the flrst derivatives. 
This fact gives rise to several technical problems. In particular, the potential of the dynamical system is no 
longer polynomial, but presents poles for some values of the variables. 

In this paper, we use the methods of the theory of dynamical systems to investigate the case of EH-GB 
theory coupled to scalars. The interest of this study relies both in the possibility of classifying EH-GB black 
holes and in the demonstration of the applicability of these methods also to higher-derivative models of the 
Lovelock type. 

The result of our investigation is that in four dimensions only asymptotically flat regular black hole 
solutions can exist, while for D > 4 we are not able to exclude the possibility of more general asymptotics. 

2. The dynamical system for 4-dimensional spacetime 

We start by considering the four-dimensional case and will then generalize to higher dimensions. Since 
in four dimensions the GB term is a total derivative in absence of scalars, the calculations are simpler, and 
the results may be qualitatively different from those obtained in higher dimensions. 

The four- dimensional action is 

I = J d'x [n - 2{d^f + ae-^'^S] , (2.1) 

where <S = TZ^'^^'^TZ^i^pcr — ^TZ^^TZfjn, -\- is the GB term, is a scalar fleld and a is a constant with 
dimension (length)^. This action also emerges as a low-energy expansion of string theory [3,5]. In this case 
the constant a is proportional to the slope parameter a' . 
We adopt the spherically symmetric ansatz [15] 

rf.S^ = -e2.^^2 ^ ;^-2g2(2C-.)^^2 ^ e2(C-.)^^2^ cj) = V - X, (2.2) 

An independent scalar charge is however present in the case of multiple scalar flelds [13]. 
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where v, (, x and A are functions of ^ and df22 is the volume element of a 2-sphere. 

Substituting the ansatz (2.2) into (2.1) and performing some integrations by parts, the action can be 
cast in the form 



xi A 



(2.3) 



As usual, the action (2.3) does not contain derivatives of A, which acts therefore as a Lagrangian multiplier 
enforcing the Hamiltonian constraint. Another relevant property of the action is that, in spite of the presence 
of the higher derivative GB term, after integration by parts it contains only first derivatives of the fields, 
although up to the fourth power, and therefore gives rise to second order field equations. 

We vary the action and then write the field equations in first order form in terms of the new variables. 



X = C, Y = v', W = x!, Z = e^, V = ^e^, 



which satisfy 



Z' = XZ, V' = WV. (2.4) 
Varying with respect to A and then choosing the gauge A = 1, one obtains the Hamiltonian constraint 



E = 2Y^ - X'^ + W'^ - 2YW + Z'^ + 8Y{Y -W) Z'^ - 3{Y - Xf 
while variation with respect to v, ( and x gives the other field equations 

2F' - W"' + 4 1 {2Y -W)Z^ -{¥ - X){4Y'^ - 2XY + XW - 3YW) 
X' -si^YiY -W){Y - X)^^ = Z'^ + 8Y{Y -W) Z'^-2{Y-Xf 



(2.5) 



Y' - W + A{Y 



Z^-{Y- X)2j ^ \ = -8Y{Y -W) Z^-{Y- Xf 



Yl' 
Yl 

z^' 



0, 



(2.6) 
(2.7) 
(2.8) 



In these variables, the problem takes the form of a five- dimensional dynamical system, subject to a constraint. 
Notice that the function E defined in (2.5) is a constant of the motion of the system, whose value is set to 
zero by the Hamiltonian constraint. 

When the coefiicient of V^/Z^ in (2.5) does not vanish, one can solve for V'^/Z'^ in terms of the other 
variables as 

(2.9) 



where 



Z^ 8Y{Y - W) [Z2 - 3(y - X)2] ' 



p2 = -X'^ + W'^ - 2YW + 



(2.10) 



and then substitute in the other field equations. In this way one obtains a four-dimensional system, with 
equations 



2Y' -W' - 



{2Y - W)Z'^ -{y - X)(4y2 _ 2XY - ZYW + XW) 



2Y{Y - W)[Z'^ - 3{Y - X)2] 



= 0, 



X' 



(Y^X) 



-3{Y -xy 



z' 



Y' -W 



Z^-{Y-Xf 



Z^-2{Y-Xf , 

z2-3(r-x)2 ' 



' Z^-{Y-X) 



2(y- w)[Z2-3(y-x)2] 



pz 



Z2 - 3(F - X)2 ' 



(2.11) 
(2.12) 
(2.13) 



z' = xz. 



(2.14) 



The substitution (2.9) allows one to avoid the singularity of the differential equations at Z = 0, but introduces 
new singularities at the poles of (2.9), that will complicate the discussion of the phase space. 
The system may be put in canonical form solving (2.11)-(2.13) for the first derivatives, as 



X' = F{X, Y, W, Z), Y' = G{X, Y, W, Z), W = H{X, Y, W, Z) . 



(2.15) 



The expressions so obtained are awkward and shall not be reported here. 

Another possible simplification of the system arises from the fact that, in analogy with the EH theory, 
the action depends on u only through its first derivatives, and therefore its variation (2.6) with respect to v 
determines a first integral, that can be written as: 



{2Y - W)Z'^ -{y - X)(4y2 _ 2XY - ZYW + XW) ^ _ 

2Y{Y-W)[Z^-3{Y-X)^] -P -a. 



(2.16) 



One may exploit this fact to reduce the problem to the study of a three-dimensional dynamical system, 
by eliminating another dynamical variable, say W, by means of (2.16). However, (2.16) gives rise to a 
complicated third order algebraic equation, that may have more than one solution for given values of X,Y 
and Z. For this reason we prefer to consider the phase space as a three-dimensional surface embedded in a 
four-dimensional space. 

In the discussion of the dynamical system, the knowledge of the exact solution for a = 0, corresponding 
toV = 0, will prove useful. In fact for a = the GB term is decoupled and hence, by the no-hair theorems, 
one should recover the usual Schwarzschild solution. In this case, the field equations reduce to 



v" = 0, 



Their integration yields 



2oe"« 
1 - e2«« ' 
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with a, b, c integration constants. Substituting in the constraint (2.5), one gets — 26^ + 26c — = 0. 
If one also requires that the solution has a regular horizon, one must impose a = b, which, combined with 
the previous constraint, yields a = b = c. Substituting these results in the metric (2.2) and defining a new 
coordinate r = 2a/(l — e^"^), one obtains the Schwarzschild metric with mass a and constant dilaton: 

ds^ = -(1 - 2a/r)df + (1 - 2a/r)-^dr'^ + r^dVt^, cj) = const. 

In the variables of the dynamical system the solution reads 

X = ocotha^, Y = a, W = a, Z 

For a — > 0, the solution reduces to flat space. 

3. The phase space of the 4-dimensional solutions 

The discussion of the phase space is complicated by the presence in the differential equations (2.11)- 
(2.13) of a denominator that diverges for F = 0, F = W, or Z"^ = 3{Y — X)^, corresponding to infinite 
or undetermined values of V. The hypersurfaces defined by these conditions are therefore singular for the 
trajectories. Although the solutions of interest do not lie on these hypersurfaces, their endpoints may. One 
must therefore be careful in the discussion of critical points. It must also be noticed that the Z <0 region 
is simply a copy of the Z > region, and hence we shall not discuss it. 

The critical points of the dynamical system are defined as the points where the derivatives of all the 
variables vanish. Setting to zero the r.h.s. of (2.12)-(2.14), one obtains the following solutions: 

I) Zo = 0, Fo = 0. 

The constraint (2.16) then implies Wq ^2Yo-a and Xq = ±^yW^~2aYoTa^. From (2.9), it follows that 
Vo = 0. This set of critical points lies on a hyperbola in the X-Y plane and on a straight line in the W-Y 
plane. 

II) Xo = 0, Fo = Wo = ±Zo. 

These points are placed on the singular hypcrsurface Y — W, and in fact V ^ oo there, except when 
Yb = Wo = Zq = 0. Therefore, only this case must be considered. It coincides with the limit a ^ of case I. 

The trajectories of the dynamical system describe the behavior of the solutions between two critical 
points. Critical points at finite distance in phase space can correspond either to singularities or to horizons 
of the physical metric. At a horizon, e^" vanishes, while the radius R = e'^~" of the 2-sphere takes a finite 
value. Instead, a critical point with = corresponds in general to a singularity. Moreover, we require that 
also the dilaton is regular at the horizon. This happens if e~^"^ = e^''^""^ is finite there. 

These conditions are equivalent to the requirement that at the critical point F^, where 

{/2 = e^"^. The critical points at finite distance are attained for ^ — > — oo, if Xq, Yq, Wq > 0, or for ^ — > +00, 
if Xo, Fo, Wo < 0, and the metric functions behave there as ~ e^^"?, U"^ ~ e^^^i, V'^ ~ e^^"^. It follows 
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a 

sinh 



that the critical points describe a regular horizon rather than a naked singularity iff Xq = Iq = Wq = a. 
It must be noticed that these points lie on the singular hypersurfaces, so one must be careful in taking the 
limits when approaching them. 

The behaviour of the solutions near the critical points can be studied by linearizing the field equations. 
Unfortunately, as already noticed, in this limit (2.9) becomes undetermined and therefore the phase space 
is effectively five-dimensional near these points. As a consequence, the linearized equations depend on how 
the limit X = W = Y is taken. A more rigorous treatment must be performed in the full five- dimensional 
phase space, and will be exposed in more detail in the appendix. 

Prom this investigation it results that for a > the critical points act as repellors, while for a < 
they attract the trajectories. In the limit a = all eigenvalues are null, and correspond in fact to constant 
solutions. 

For the investigation of the asymptotic behaviorof the solutions it is useful to study the critical points 
at infinity [15]. In fact, spatial infinity is achieved for R = e^~" — > oo, i.e. e^ = Z — > oo. This limit, which 
in general is reached for a finite value of ^, corresponds to the surface at infinity of phase space, where 
the fields diverge. 

Phase space at infinity may be studied by defining new variables 

I Y W Z P 

M = ^, y = ^, "'^x' ^^x' ^"x" ^ ' 

The field equations at infinity are then obtained for u ^ and may be written as 

u = -fu, y = 9- fy, w = h- fw, i = (1 - f)z, (3.2) 

where a dot denotes ud/dS, and / = u^F{l, y, w, z), g = m^G(1, y, w, z), h = u^H{l, y, w, z). The constraint 
(2.16) reduces to 

^ ""^ 2y{y - w)[z^ - 3{y - m ^ " ^^'^^ 

The critical points at infinity correspond to the vanishing of the derivatives in (3.2) for u — > 0. This 
occurs for 

I) yo = Wo — 0, zo — ±1, i.e. Xq = ±oo, Yq — Wq — 0, Zq — ±Xo. In this case vq is undetermined. 
These are also the endpoints of the trajectories lying in the V = hypersurface, corresponding to the exact 
solutions discussed at the end of the previous section. 

II) zo = 0, po = 0, that in view of (3.3) yields yo = "^f i ''^o = V^. These points correspond to 
Xq = ±00, Fo = ^7 Wq = y/2Xo, Zo = 0. It follows that vq = 0- These are the endpoints of the 
hyperboloid Z^ = P"^ = 0. 

Also in this case, phase space is five-dimensional near the critical points, and its linearization will be 
studied in more detail in the appendix. It will result that for Xq > 0, points I attract all the trajectories. 



while points II repel the trajectories at infinity, and are centers for those coming from finite distance. For 
Xo < repellors become attractors and viceversa. 

Prom the location of the critical points at infinity, one can deduce the asymptotic behaviour of the 
solutions [15]. In case I, C ~* ^o, and 

^ — ^o|~"^) ^ ^ const, ^ const, 

and therefore, defining r ^ |^ — Co|~^ and recaUing (2.2), one obtains after some calculations the asymptotic 
behavior 

ds^ ~ -dt^ +rfr^ +r^(if^2i ~ const. 
Hence these points correspond to asymptotically flat solutions with asymptotically constant scalar field. 
At points 11, ^ ±00, depending on Xq being ±oo, and 



and therefore, defining r e^^^. 

Summarizing, from the study of phase space emerges that all regular spherically symmetric black hole 
solutions start from the point Xq = lb = M^o = a > 0, = and end at the critical point 1 (with X > 0) 
at infinity (or start at the critical point I at infinity with X < and end at Zq = 0, Xq = Iq = ^'^o = a < 0). 
It follows that all solutions are asymptotically flat, with asymptotically constant scalar field. In particular, 
no solutions with anti-de Sitter asymptotics exist, in accordance with the results of [14]. 



4. The dynamical system in (m + 2)-dimensional spacetime 

In this section we generalize the previous results to the case of m + 2 dimensions. Since the procedure 
is identical to the case m = 2, we only give the main formulas. 
The (m + 2)-dimensional action can be written as 

7 = y d^^+^^a: [7^ - m(m - \){d(l>f + ae-'^'^S] . (4.1) 

We adopt the spherically symmetric ansatz 

rfs2^_g2(,n-l).^^2^;^-2g2(mC-.)^2^g2(C-.)^j^2^^ = - (4.2) 

Substituting in the action and after integration by parts, one has 

p2(— DC + (m - 2) (m - 3)a e2x+2(m-2)c 



/= m(m-l) / ^(^+2)2; \ - 



-A 



2!/'2 - + x" - 2!/'x' 



-2QAe2x-2C 



(m^ - m + 2)!/'2 



- (m - 2)(m - 3)C" - ^v'y^ - 4(m - 2)C'x' 



+ ^ e2'^-2"C(C' - v'f 3(m2 + 3m - 2)1/'^ _ (^ - 2)(m - 3)C'^ 



2(m - 2)(m - 3)CV - 8(2m - l)i/'x' 



2)CY 



(4.3) 



As for m = 2, the higher- dimensional action (4.3) has the property of containing only first derivatives of v, Q 
and X and no derivative of A. However, more terms appear in the action, because the GB term is no longer 
a total derivative in absence of scalar fields. 

We vary the action and write the field equations in terms of the new variables. 



which satisfy 



Z' = {m- 1)XZ, V' = WV. (4.4) 
Varying with respect to A and then putting A = 1, one obtains the Hamiltonian constraint 

p2 + |(m - 2)(m - 3)^4 + 2 f (m^ - m + 2)Y'^ - (m - 2)(m - 3)^^ - 4YW 



(4.5) 



4(m - 2)XW Z^-{X- Yf 3(m^ + 3m - 2)Y^ - (m - 2)(m - 3)^^ 

y2 



2(m - 2)(m - S)XY - 8(2m - 1)YW - 8(to - 2)XW 



^2m/(m-l) 



= 0, 



where is given by (2.10). 

Variation with respect to u, C, X) respectively, gives the other field equations 

2Y'-W' + 2l^ ((m^ -m + 2)Y- 2W^ Z'^ + {X - Y) ((m^ + 3m - 2)Y'^ 

-{m'-m + 2)XY - 2(2m - l)YW + 2Xw)^ ^2^^] = 0' 

X' + 2| (m - 2) ((m - 3)X + 2W^ Z'^ + ^{X -Y)(^ 2m(m +l)Y^-{m- 2)(m - 3)^^ 



- (m - 2)(m - 3)XY - QmYW - 6(m - 2)Xw) ^ ^2^(^-1) 
(m - 1)^2 + (m - 2f{m - 3)Z^ + 2((m^ - m + 2)Y'^ - (m - 2)(m - 3)^^ - 4YW 

- 4(m - 2)XW^ Z'^-j{X- Yf ( 3(m2 + 3m - 2)Y^ - (m - 2)(m - 3)^^ 

- 2{m - 2)(m - 3)XY - 8(2m - 1)YW - 8(m - 2)XW] 



^2m/(m— 1) ■ 



Y' -W' + 4 



( r + (m - 2)X) - i(X - Yf ((2m - 1)F + (m - 2)X 
(m - 2)(m - 3)Z^ - 2(^{m'^ - m + 2)y2 - (m - 2)(m - 3)^2 - AYW 

- 4(m - 2)XW) Z2 + i (X - y)2 ( 3(m2 + 3to - 2)y2 _ _ 2) (m - 3)X 

- 2(m - 2)(m - 3)XY - 8(2m - 1)YW - 8(m - 2)XW] 



2'2m/(m-l) 



(4.6) 



(4.7) 



^2m/{m-l) ' 



(4.6 



Of course for m = 2 the equations reduce to those of the previous section. Again, (4.6) gives rise to a 
first integral: 



2Y -W + 2 



{m^ -m + 2)Y - 2W Z'^ + {X - Y) [m^ + 3m - 2)Y^ 



-[m''-m + 2)XY + 2XW - 2(2m - l)Yw] | ^^^^^ = a. (4.9) 

The factor y2^^2m/(m-i) ^j^g equations (4.6)-(4.8) can be eliminated by use of (4.5). In this way 
one obtains rational expressions in the variables X, Y, W and Z. The equations can then be cast in the 
canonical form 

X' = F{X,Y,W,Z), Y' = G{X,Y,W,Z), W = H{X,Y,W, Z). (4.10) 

As before, the first integral (4.9) may be used to eliminate from the system the variable W and reduce 
the problem to a three-dimensional dynamical system, but we prefer to consider the phase space as a 
three-dimensional surface embedded in four-dimensional space. Also, one must be careful with the possible 
singularities introduced by solving (4.5) with respect to \/2^^2m(m-i)^ 

The calculations were performed with MATHEMATICA. Unfortunately, the expressions of the functions 
F, G and H contain hundreds of terms and remain unmanageable also using computer algebra: our analysis 
is therefore necessarily limited by this problem. 

The critical points of the dynamical systems at finite distance are 

I) Zq = 0, Pq = 0, and hence Vq = 0. The constraint (4.9) then implies Wq = 2Yq — a and Xq = 
±^/2YQ — 2aYo + a^. These points correspond to a regular horizon only \i Xq = Yq = Wq = a. 

II) A generalization of the point II in m = 2, which is defined by Xq = 0, Zq = 2YoWo — 2Yq — Wq, with 
Yo and Wq related by the equation (m - 2) (to - 3)W^ - 4(to - 1) (to - 4)Wf^Yo + 2{3m'^ - 19m + 14)W^Y^ 
-|(m^ - 23m + 16)VKoi^o^ + {m^ - 13m + 14)1^0^ = 0. As for m = 2, these points correspond to a finite 
value of Vo only if Yq — Wq — 0, where they reduce to a special case of I. 

The phase space at infinity may be studied like for m = 2 by defining the new variables u, y, w and z 
and taking the limit u — > 0. The field equations at infinity can then be written as 

u = -fu, y = g-fy, w = h-fw, z = {m-1- f)z, (4.11) 

with the notation introduced in the previous section. The critical points correspond to the vanishing of the 
derivatives. This happens for 

I) Wq = yo = 0, zq = ±1, with Vo undetermined. 

II) zo=Po = 0, and hence yo = = \/2, vq = 0. 

Due to the size of the equations involved, we cannot however exclude the presence of further critical points, 
although this seems unlikely. 
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For this reason, we do not discuss in detail the linearization of the field equations near the critical points, 
although the results look similar to the four- dimensional ones. Also the asymptotic behaviornear the critical 
points is analogous to m = 2. In case I, ^ — > Coi and 



e? ~ 1^ - ^o\-^/^"'-^\ ~ const, ~ const, 
and hence, defining r |^ — ^o\~^^^"^~^^ and recalling (4.2), 

ds'^ ~ -dt^ + dr^ + r'^dn'^, cj) ~ const. 



At points II, ^ — > ±oo, and 



and defining r ~ e=F['"+('"-2)/^)l«, 

(^^2 ^2(m-l)/[(m-2)+V2m]^^2 _^ ^-2(m-l)/[(m-2)+v^m]^^2 ^-2(1- V2)/[(m-2)+N/2m]^Q2^^ 

g-20 ^ ^2/[(m-2)+V2m]_ 

We can conclude that regular asymptotically flat black hole solutions exist also in D > 4, but we cannot 
completely exclude the presence of further solutions with different asymptotics. Asymptotically (anti)-de 
Sitter solutions can nevertheless be ruled out by explicit check. 

5. Conclusions 

We have studied the dynamical systems generated by the equations that determine the spherically 
symmetric solutions of EH-GB dilatonic theories and shown that in four dimensions the only regular black 
hole solutions are asymptotically fiat. In Z) > 4 a complete answer was not possible due to the complexity 
of the system, but it seems plausible that the same property holds. 

In spite of the diflaculty of studying the case > 4, we find however important the demonstration that 
methods from the theory of dynamical systems can be used also in the case of higher-derivative models of the 
Lovelock type. In particular, these methods may be used to study more involved problems, as for example 
black hole solutions arising from the compactification of higher dimensional EH-GB models, or cosmological 
models. 

Appendix 

In this appendix we give some details of a more rigorous treatment of phase space in the four- dimensional 
case. This can be achieved if, instead of confining to the physical hypersurface E = 0, one considers the full 
phase space. In this way, the only singularity is at Z = 0, and it is easier to control the behaviour of the 
trajectories near the critical points. 

Since we are not aware of any treatment of singular dynamical systems in the literature, the only 
possibility to proceed rigorously would be to define new variables such that the system becomes regular, 
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for example introducing V = V/Z"^. These variables however are not suitable for discussing the physics, 
since V would diverge in physically interesting situations, as for example at the horizon of Schwarzschild-like 
solutions. For this reason we shall attempt to maintain the old variables, taking the limit Z ^ in the 
correct way, when ambiguities may occur. 

Our starting point are equations (2.4)-(2.8). Since E in (2.5) is a first integral of the dynamical system, 
all the trajectories are confined to surfaces of constant E, but only those with E = are physical. We shall 
therefore consider only the critical points lying on this surface. Equations (2.6)-(2.8) can be explicited for 
X', Y' and W, as in (2.16): 

X' = F{X,Y,W,V,Z), ¥' = G{X,Y,W,V,Z), W' = H{X,Y,W,V, Z), (Al) 

where F, G and H are rational functions of their variables. 

Let us consider the critical points at finite distance first. These are placed at = Vq = 0. To simplify 
the discussion, we shall only consider the points corresponding to regular horizons, where Xq = Yq = Wq = a, 
and 

Z ~ y ~ e°«, {A.2) 

with ^ — > ±oo, depending on the sign of a, as explained in sect. 3. In view of (A.2), we linearize the equations 
(2.4) and (2.6)-(2.8) around the point Xq = Yq = Wq = a, Vq = kZo — > 0, for finite k, obtaining 

SY' =6X' = 6W' =4a6Y, 6Z' = adZ, SV' = a6V, (A3) 

with eigenvalues (2), a (2), 4a, independent of k. The zero eigenvalues are due to the presence of two 
constants of motion. It results that the critical point is an attractor for a < and a repellor for a > 0. 
As discussed in sect. 3, the phase space at infinity can be studied by defining new variables 

1 Y W Z V 

u = ^, y=^, w = -, z = -, v = -. 

The field equations at infinity (w — > 0) read 

u = -fu, y = g-fy, w = h-fw, z = {i-f)z, v = {w-f)v, 

and the critical points are of course those obtained in section 3: 

I) yQ = wo = 0, zq = ±1, any vq, i.e. Xq = ±oo, Yq = Wo = 0, Zq = ±Xq, any Vq. 

II) Zq = 0, Vq = 0; for the physical solutions yo = wq = ±\/2. In terms of the original variables, 
Xo = ±00, Yo = ^, Wo = V2X0, Zo = Vo = 0. 

The linearization around points I poses no problem, yielding 

dii = —du, Sy = —Sy, dw = —6w, 6v = —6v, Sz = —26z. (^-4) 

The eigenvalues are therefore all negative and for Xo=oo (Xo=— oo) these points attract (repel) the trajec- 
tories both from finite distance and from infinity (remember that ii = ud/d^). 
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At points II, for ^ — > ±00, Z ^ e^^, V ~ e^*^^, and hence vjz — > 0. Therefore, one must take the limit 
u — > first, and then 2; — > 0, getting 

Sii = 6y = Sw = 0, 6v = V2 6v, Sz = Sz. (^-5) 

The eigenvalues are either null or positive. These points are therefore centers for trajectories coming from 
finite distance, and repellors (actractors) for those with Xo=oo (Xo=— 00). 

The resulting phase space portrait is rather similar to that of the EH-scalar system, except for the 
location of the critical points II at infinity, that in the EH case are placed at yo = wo = 0. In particular, all 
the trajectories corresponding to solutions with regular horizons start from the critical point a,t Z = V = 0, 
X = y = W^ = a>Ofora parameter a, ending a,t Z = X = 00, Y = W = 0, and hence are asymptotically 
fiat. These solutions also possess ^-reversed copies that start at Z = X = —00, Y = W = 0, and terminate 
a,tZ = V = 0, X = Y = W = a<0, with identical properties. 
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